Spectral statistics of disordered metals in the presence of several Aharonov-Bohm 

fluxes 
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The form factor for spectral correlations in a diffusive metal is calculated in the presence of several 
Aharonov-Bohm fluxes. When the fluxes <f) are equal, the correlations are universal functions of 
ng 2 (f> where g is the dimensionless conductance and n is the number of applied fluxes. This explains 
recent flux dependence of the correlations found numerically at the metal-insulator transition. 
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In the recent past, there has been several studies on 
the effect of a Aharonov-Bohm flux on the spectral prop- 
erties of disordered metals. These works were mainly 
stimulated by the observation of persistent currents in 
mesoscopic rings . In addition, the flux can be used as 
a tool to probe the evolution of the statistical properties 
of spectra between different symmetry classes 
When no flux is applied, time reversal symmetry is pre- 
served, the low energy spectral correlations are those of 
an ensemble of real random matrices [^|,0J^1 . This ensem- 
ble is invariant under orthogonal transformations and is 
called the Gaussian Orthogonal Ensemble (GOE). In a 
finite flux, time reversal symmetry is broken and the rel- 
evant ensemble of random matrices is the Gaussian Uni- 
tary Ensemble (GUE). One or several magnetic fluxes 
drive the transition between the two ensembles. A mag- 
netic flux is also an example of external parameter which 
can probe the parametric correlations, i.e. correlations 
of energy levels between different values of a parameter 

Quite recently, a Aharonov -Bohm flux has been ap- 
plied in numerical experiments to study the transition 
between universality classes at the metal-insulator tran- 
sition 1 1 1 In this work, several fluxes were applied in 
perpendicular directions and the authors noticed that the 
cross over depends on the number of applied fluxes. 

The aim of this short note is to describe analytically 
the effect of several fluxes. Namely, I extend the semi- 
classical calculation of the spectral correlations in metals, 
and in particular of the form factor, to the case where 
several fluxes are applied. Consider a three-dimensional 
disordered metal in the diffusive regime, i.e. when all the 
dimensions L Xl L y and L z are larger than the mean free 
path l e , so that the electronic motion is diffusive in the 
three directions. Periodic boundary conditions are im- 
posed in the three directions but the hypertorus can be 
pierced by magnetic fluxes <p x ,<f)y, 4>z in three perpendic- 
ular directions. It is known that the effect of a flux line 
is simply to change the boundary conditions so that: 



^(x + L x ,y,z) = e 12 ^ 
i;(x,y + L y ,z) = e' i27r ^ rp(x, y, z) 
ip(x,y,z + L 

<J)i/(j)Q, with i 



ip(x,y,z) 



ij)(x,y,z) 



(1) 
(2) 
(3) 



x, y, z. 4>q is the flux quan- 



where ipi 
turn. 

To describe the spectral correlations, I use the semi- 
classical result found by Argaman et al. Jll[ to describe 
the two-point correlation of the density of states K (e) = 
(p(E)p(E + e)) - (p(E)) 2 . These authors have related the 
form factor K (t), the Fourier transform of this two-point 
correlation function (h = 1 throughout the paper): 

K(t) = -!- / K(e) exp{iet)dt (4) 
2vr J 

to the return probability P(t) for a diffusive particle. In 
the presence of magnetic fluxes ifi, this relation is: 

K(t, tpi) — tP(t, <pi)/ (47r 2 ) . (5) 

When the phase of electrons is coherent, the return 
probability contains two terms, the classical return prob- 
ability P c i and an interference term, Pi n t, which results 
from phase coherence between a diffusive path and the 
same path obtained by time reversal symmetry. The lat- 
ter term depends on the magnetic fluxes. It is the solu- 
tion of the diffusion equation and is given by: 



Pint(t,<Pi) 



{qx,q v ,qz} 



D{ql+ql+ql)t 



(6) 



where the wave vectors of the diffusion modes are deter- 
mined by the boundary conditions: qi — 2(71; + 2(pi)ir/Li 
and Ui = 0, ±1, ±2, • • -. I? is the diffusion coefficient. For 
our purpose, it is more convenient to rewrite this inter- 
ference term as a sum over winding numbers around each 
flux tube: 



Pint (t, <fi 



L x L y L z 



E 



{rrii}=- 
X C0S(47T niilfii) 



(7) 
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Li the dimension of the system in the direction i. The 
sums are over i = x,y, z. This expression can be de- 
duced from eq. ^ by Poisson summation. Quite simply, 
this result expresses that the mf harmonic of the flux 
dependent part of the return probability depends on the 
probability to return to the origin after rm windings on 
the hypertorus. This is a simple generalization of the one 
flux case Q. Thus the spectral correlations, in particu- 
lar the form factor K(t), are related to the way the phase 
is accumulated along paths of diffusive trajectories after 
enclosing the different fluxes. 

In the limit of small fluxes ipi <C 1, the flux periodicity 
can be neglected and the above sums can be replaced by 
gaussian integrals. One gets: 

Pint(t, <fi) = CXp[-f67T 2 ^ E cifit] 

i 

where E c i — TiD/L 2 is the Thouless energy in the direc- 
tion i. The form factor is thus given by: 

The first term of this expression comes from the classical 
return probability which is equal to the interference term 
when the fluxes are zero. 

When all the dimensions and applied fluxes are identi- 
cal, 4>i — (j), E ci — E c , a situation considered in ref. ]To| |, 
the form factor depends on the single combination of pa- 
rameters nE c ip 2 t. This implies that the flux dependence 
of the two-point correlation function is solely governed by 
ngip 2 where g = E c / A is the dimensionless conductance, 
n is the number of applied fluxes and A is the mean level 
spacing. As expected, the transition between orthogonal 
and unitary symmetries is thus faster when three fluxes 
are applied. 

Our calculation has been done in the semiclassical 
regime where the time is smaller than the Heisenberg 
time tji = h/A. The form factor for the orthogonal- 
unitary transition is completely known from Random 
Matrix theory |§ or from the non-linear a model ||. 
More generally, one expects all correlations functions to 
be universal functions of the combination ngtp 2 . 

Let us now come to the spectral correlation at the 
metal-insulator transition. At the transition, the con- 
ductance g is scale independent so that the transition 
between orthogonal and unitary symmetry classes is ex- 
pected to be size independent. This scale invariance has 
been found numerically recently by Batsch et al. JTo[ ] . In 
their work, they study the evolution of the distribution 
P{s) of spacings between nearest energy levels. They in- 
troduce the parameter T max (ip) = (P m (<p) ~ P r ° n )/(P™ - 
Pf^), where P m is the maximum of the spacing distribu- 
tion for the orthogonal (Pm), unitary (P^) and interme- 
diate (P m {tp)) cases. To describe the cross-over between 
orthogonal and unitary symmetries, they present differ- 
ent plots of the variation T max (tp), for the three cases 



where one, two or three identical fluxes are applied. They 
find that each flux dependence is scale invariant at the 
transition, in accordance with the scale invariance of the 
conductance. However, they obtain three different vari- 
ations (fig. 3 in ref. E(J), depending on the number of 
applied fluxes, a result that they qualify as "unexpected 
behavior" . 

As for K(e, tp), one expects the quantity T max (ip) to be 
a universal function of ngtp 2 in the metallic regime. One 
may now wonder if this scaling also holds at the metal- 
insulator transition. If it is so, one expects that the three 
curves obtained in fig. 3 of ref. [[lOj, for one, two and three 
fluxes, merge into a single universal curve by rescaling 
ip — > ip/y/n. This is indeed the case as it is shown on 
fig. I. I have used the data of the fig. 3 in ref.[ p0|]. The 
flux scales have been renormalized as <p — > p\i for one 
flux, ip — > ipy/2/y/3 for two fluxes, so that all the curves 
now coincide with the three-fluxes curve. This shows that 
indeed T max (ip) is a unique function of rup 2 . 

It may appear surprising that the scaling found in the 
diffusive regime still applies at the Anderson transition. 
At the transition, the correlation functions have indeed 
a different structure. In particular, it has been found re- 
cently that the form factor is a function of ip 2 t 1+r, ^ d where 
r\ is an exponent related to the multifractal structure of 
the wave function at the transition Jl3| |. However, several 
fluxes still contribute additively to the total accumulated 
phase so that the orthogonal-unitary cross-over is still a 
unique function of ng*(p 2 at the transition, where g* is 
scale invariant. 

In conclusion, we have calculated the spectral form 
factor in the presence of several Aharonov-Bohm fluxes. 
When these flux are equal, the Orthogonal-Unitary tran- 
sition is driven by the combination of parameters ngip 2 . 
This result still holds at the transition and explains re- 
cent numerical results. 

I thank M. Pascaud and P. Walker for useful com- 
ments. 
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FIG. 1. The universal function T max (ip) after the rescaling 
proposed in the text. The two small dots, which correspond 
to the case of one flux are outside the universal curve, because 
the transition cannot be completed before the flux periodicity. 
For the legend of the symbols, see ref. [juj. The continuous 
line is a fit to the data. 
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